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Abstract. The dynamical density functional theory of Marconi and Tarazona [J. 
Chem. Phys., 110, 8032 (1999)], a theory for the non-equilibrium dynamics of the one- 
body density profile of a colloidal fluid, is applied to a binary fluid mixture of repulsive 
Gaussian particles confined in a spherical cavity of variable size. For this model fluid 
there exists an extremely simple Helmholtz free energy functional that provides a 
remarkably accurate description of the equilibrium fluid properties. We therefore use 
this functional to test the assumptions implicit in the dynamical density functional 
theory, rather than any approximations involved in constructing the free energy 
functional. We find very good agreement between the theory and Brownian dynamics 
simulations, focusing on cases where the confined fluid exhibits phase separation in the 
cavity. We also present an instructive derivation of the Smoluchowski equation (from 
which one is able to derive the dynamical density functional theory) starting from the 
Liouville equation - a fully microscopic treatment of the colloid and solvent particles. 
This 'coarse-graining' is, of course, not exact and thus the derivation demonstrates the 
physical assumptions implicit in the Smoluchowski equation and therefore also in the 
dynamical density functional theory. 



1. Introduction 

For a multi-component fluid composed of particles with a large size difference between 
the different components, a theoretical description of the (inhomogeneous) fluid 
dynamics is challenging due to the varying equilibration timescales for the different 
species of particles in the fluid. Colloidal suspensions are a particular example of 
such complex fluids. Typically, the colloids are 10-1000nm in diameter, whereas the 
solvent molecules are, in the case of water, ~ 0.2nm. The situations in which a theory 
could be used to describe the motion of inhomogeneous (confined) colloidal particles 
are many. For example, in colloidal physics the use of optical tweezers to confine and 
then move individual of groups of colloids is common pL . Simulating such a fluid is also 
computationally expensive, because of the huge numbers of solvent molecules required 
when simulating even a limited number of colloids. A common approach to such systems 
is to integrate out the degrees of freedom of the smaller particles in order to render a 
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description of the system based only on the (much smaller) phase space corresponding 
to the degrees of freedom of the big colloid particles [2 . 

In the case of colloidal fluids in thermal (static) equilibrium, this 'coarse-grained' 
approach is, at least formally, well understood and one is able to formulate an effective 
Hamiltonian for the fluid involving only the phase space coordinates of the big particles. 
Of course, the effective potential between the big particles is, in general, many-body 
in character and dependent on the density of both the solvent and colloid particles 
[21 Ei|- For example, the depletion potential, or effective solvent mediated interaction 
potential, between big hard-sphere particles dissolved in a solvent of small hard spheres, 
is oscillatory for large separations 0] and is therefore very different from the 'bare' 
potentials one typically encounters between the particles of a molecular fluid. However, 
for the non-equilibrium dynamics of colloidal fluids, coarse-graining is carried out on a 
more ad-hoc basis. A common approach is to use the effective solvent mediated potential 
obtained from the equilibrium theory together with stochastic equations of motion - i.e. 
the force on the colloids due to collisions with the solvent particles is modelled as a 
random white noise term and a frictional (Stokes) one-body drag force. The colloids 
are thus modelled as Brownian particles. From the (Langevin) equations of motion for 
the Brownian particles one can obtain the Fokker-Planck (Smoluchowski) equation [3] 
for the colloid probability distribution function in phase space, and thus one is able to 
determine non-equilibrium dynamic properties of the colloidal fluid. In going to such 
a stochastic description of the dynamics one inevitably neglects hydrodynamic effects. 
We justify the above stochastic approach in Sec. |H 

In this paper we consider cases of a model colloidal fluid confined in time dependent 
one-body external potentials. We focus on the time dependence of the fluid one- 
body density profile, using the dynamical density functional theory (DDFT) derived 
recently by Marconi and Tarazona |H1 Ej §■ In a recent paper i9| (see also Ref. 
[Til]), it was demonstrated that this DDFT could be derived from the Smoluchowski 
equation, by making the assumption that the form of correlations between particles in 
an inhomogeneous fluid out of equilibrium, is the same as in an equilibrium fluid with 
the same one body density profile [HJ E]. We apply the DDFT to the dynamics of 
a particular model colloidal fluid, the binary Gaussian core model (GCM), making 
a comparison between the results from the DDFT and Brownian dynamics (BD) 
simulation results. Dzubiella and Likos [TT] made a similar comparison for the one 
component GCM, and found good agreement between theory and simulations. The 
GCM is a simple model for the effective potential between the centers of mass of polymers 
or dendrimers in solution [3 US HH1 EU C3 El EZl HE] • The reason Dzubiella and Likos 
chose this particular model fluid is that there exists a very simple, yet remarkably 
accurate (as we demonstrate below for mixtures) Helmholtz free energy functional for 
describing the equilibrium properties of the GCM HH1 H3 EU E21 EH1 E3 and thus 
by using this functional combined with the DDFT they were able to test some of the 

§ The key DDFT equation (Eq. HI 60 ) was proposed originally, without derivation, by Evans (HJ- 
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approximations inherent in the DDFT, rather than the approximations involved in the 
free energy functional. The DDFT has also been successfully applied to situations with 
steady currents j2SJEH] and very recently, Rex et al have applied the DDFT to the GCM 
fluid in shear and travelling wave potentials j2Zj- For these cases they also find good 
agreement between theory and simulations. One particular question we wish to address 
here is whether the DDFT can incorporate dynamical effects that are specific to colloidal 
fluids composed of two different species of colloids and, in particular, dynamical effects 
arising when there is phase separation and wetting phenomena in a confined colloidal 
fluid. 

The paper is arranged as follows: In Sec. El we introduce and briefly describe the 
DDFT. The results from the application of the DDFT are presented in Sec. H3 In Sec. 
0] we give an alternative derivation of the Smoluchowski equation. Starting from the 
Liouville equations for a mixture of colloid and solvent particles, we first obtain the 
Kramers equation and then from this the Smoluchowski equation. These equations can 
also be derived as the (generalised) Fokker-Planck equations for a system of Brownian 
particles with stochastic equations of motion. The present derivation elucidates some of 
the physical assumptions concerning the fluid dynamics that are implicit in the Kramers 
and the Smoluchowski equations. We believe this derivation sheds light on the status of 
the DDFT and the approximations inherent in this theory. Finally, in Sec. Owe discuss 
our results and draw some conclusions. 

2. The DDFT 

Before introducing the DDFT, we recall briefly some of the basic results from equilibrium 
density functional theory (EDFT) jHJ EE| that are required for our discussion of the 
DDFT. The key quantity in EDFT is the Helmholtz free energy functional: 



which is a functional of p(r), the fluid one-body density. For simplicity we consider a 
one-component fluid. u ex t(r) is a one body external potential and 



is the (exact) ideal gas free energy; A is the thermal de Broglie wavelength of the 
particles. F ex [p{v)] is the excess (over ideal) contribution to the free energy due to 
interactions between the particles [HI EE] ■ In EDFT one considers the following grand 
potential functional 



where /i is the chemical potential. The equilibrium density distribution is given by the 
minimisation condition [5] I28j: 




(1) 




(2) 




(3) 



o. 



(4) 
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- rr r n - (6) 



Thus, using Eqs. (HJ) - (@J), one obtains the following Euler-Lagrange equation for the 
equilibrium fluid density profile: 

= k B T In A 3 p(r) - k B Tc {1) (r) + u ext {v) - (5) 
where the one body direct correlation function: 
JD M _ ^F ex [p(r )} 
5p(r) 

(3 = 1/ksT is the inverse temperature and — fc^Tc' 1 ^) is an effective one body potential 
that incorporates the effect of the interparticle interactions in the fluid. The central 
task in EDFT is to find a suitable approximation for F ex [p(r)] for the fluid of interest. 
F ex [p(r)] is, in general, an unknown quantity. A second functional derivative of F ex [p(r)] 
gives the inhomogeneous (Ornstein-Zernike) pair direct correlation function: 

Often, approximations for F ex [p(r)} are constructed by requiring that the approximate 
excess Helmholtz free energy functional generate, via Eq. (|7jl. reliable results for c^ 2 \r), 
the bulk pair direct correlation function, about which much is known from the theory of 
bulk equilibrium fluids |2B1E2|. 

We now move on to consider a non-equilibrium fluid of A" Brownian colloid particles. 
We denote the position of the i th colloid particle by r^, and the set of all position 
coordinates by = {ri, r N }. The total colloid potential energy Un(t n , t) is assumed 
to be of the following form: 

N I N 1 N 

U N (r N ,t) = ^2uext{Ti,t) + ^^^rj-) + - $^5Z n3 ( ri ' r J' rfc ) + ( 8 ) 

i=l j^i i=l k^j^i j^i i=l 

which is made up of a one-body term, i.e. a time- dependent one-body external potential 
u ext (ri,t) acts on each fluid particle; a two-body term, which is a sum of (time 
independent) pair potentials W2(rj,rj); a three body term, given by a sum of three- 
body potentials u 3 (ri, r J; r k ) and higher body terms, each given by a sum of multi-body 
potentials. Since we do not include explicitly the solvent particles, this potential must 
include contributions from solvent mediated interactions between the colloids as well 
as the direct (bare) interactions between the colloids - we will return to this issue in 
Sec. |H|. On the Brownian time scale the equations of motion of the colloids obey the 
following stochastic (Langevin) equations of motion [30] : 

at OYi 

|| For the purposes of calculating equilibrium fluid properties, all solvent effects are included in the 
effective potential Un(? n , t). However, for dynamics, the solvent friction has to be explicitly taken into 
account. 
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where T 1 is a friction constant characterising the one-body drag of the solvent on the 
colloidal particles and G«(t) = £f (t)) is a white noise term with the property 

(cm = o, 

{i:m(t'))=2k B T5^5(t-t'). (10) 

The stochastic equations of motion (jHJ), of course, neglect hydrodynamic interactions 
between the colloids. 

We can define a probability density function p( N >(r N ,t) for the N colloids. The 
time evolution of p( N \r N ,t) is described by the Smoluchowski equation [3]: 



dt *r~i dvi V dri dri 

i=i x 

This equation is the (generalised) Fokker-Planck equation for the Langevin equations 
(jHJ). The Smoluchowski equation is generally presented from this stochastic viewpoint. 
However, as we show below in Sec. IH one can argue for its use as an approximation to 
the exact Liouville equations: By going to the Smoluchowski equation the description 
of the fluid is reduced to one based solely on the position coordinates of the colloids, 
rather than utilising the full set of phase space coordinates for the colloid and solvent 
particles. However, in practice a further reduction is required in order to be able to 
determine explicitly the dynamics of the colloids. In particular, we focus on the colloid 
one-body density piT!HI23j: 

p(r 1 ,t)=N J dr 2 ...J dr N pW(r N ,t). (12) 

Similarly, the n-particle density is 

P W (rV) = ^^yy J dv n+l ...j dr N P^(v N ,t). (13) 

We now follow closely the derivation in Ref. [HJ: Using Eqs. (jHJ), (|T2|) and (JT3J), we find 
that on integrating Eq. (fTTj). one obtains [§1 13T]: 



^ dpjr^t) d 2 p(r u t) d 

ot orf or 1 



du ext (r 1 ,t) 



dr 1 



- ■ J dr 2 p [ ; (r 1 ,r 2 ,t) 



dv\ J dr 1 
+^ • /dr 2 J d^n, r 2 , r 3 , t) ^(w 3 ) + (u) 

At equilibrium, when dp(r,t)/dt = 0, this equation is just the gradient of the (exact) 
first equation of the YBG hierarchy |29j . 

We are now in a position to make contact with EDFT. For an equilibrium fluid 
the gradient of — ksTc^ir) (see Eq. ©), a one body force due to interactions between 
particles in the fluid, is given by the following sum- rule [HI El E2j : 

- k B T p{ll ) a -^ll = ±J dr,../ d r „p<»)( r «)^p. (15) 
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If we assume that we can use (J 15)) for a non-equilibrium fluid we are assuming that 
the effective one body force on a particle in the fluid due to interactions with the other 
particles is the same as that in the equilibrium fluid with the same one body density 
profile. We make this assumption and using Eqs. (|T5*|) and (JTJ) together with Eq. (fTlj) 
we obtain the key DDFT equation [SI 13 Ej: 



dp{r,t) _ T d_ _ 



(16) 



dt dr 

where F[p(r, t)] is given by Eq. (JTJ) with p(r) replaced by the time dependent one-body 
density p(r,t). Before commenting on the status of Eq. (|16p. we note that the above 
arguments can easily be generalised to the case where there are several different species 
of colloids. If there are Q different species of colloids, with N q colloids of species q, 
such that the total number of colloidal particles is N = Yuq=i ^<?> anc ^ the total colloid 
effective potential energy is (c.f. Eq. (jHJ)): 

Q N Q 1 Q N Q 

U N (r N ,t) =J2J2 u Ur g ,ut) + EE M " 3 '(^-V,i) + - (17) 

9=1 i=l 9:9 /= l j¥=i *=1 

In this case, the Smoluchowski equation (|TT|) simply becomes 

where T" 1 is the friction constant for the q th species of colloid particle. One can then 
integrate Eq. (fTH|) to obtain a DDFT for the one-body density profiles. The one-body 
density of species q is: 

N q 

p q (v qA ,t) = N q ]J dr q , J] / dv^'P^\v N ,t), (19) 

i=2 J q'+q J 

i.e., p q (r qi x,t) is given by the integral over with respect to all the colloid position 
coordinates other than those of the i = 1 colloid of species q. The multi-component 
generalisation of the DDFT equation (jlfi)) then becomes: 



dp q (r,t) = d_ 
dt q dv ' 



dr V Sp q (r,t) 



(20) 



where the Helmholtz free energy functional for the multi-component colloidal fluid is 
(c.f. Eq. ((H)): 

Q , 

F[{p q {r,t)}] = J2 k BT / drp q (r,t)[\n(A qPq (r,t))-l}+F ex [{p q (r,t)}} 

q=l J 



Q r 

+ Yl / dru ^* 

9=1 J 



v,t)p q (v,t). (21) 



For an inhomogeneous equilibrium fluid Eqs. (jHJ) and (jlj) imply that the chemical 
potential p = 5F/5p(r) is a constant throughout the body of the fluid. Similarly, for 
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an inhomogeneous equilibrium multi-component fluid, the chemical potentials for the 
different species, 



take a constant value throughout the fluid. The DDFT (}20|) is equivalent to assuming 
that in the non-equilibrium fluid, this is not the case, and that the gradients of the 
chemical potentials are the thermodynamic forces driving particle currents j g of each 
species IB] : 



one obtains the DDFT Eq. (|20p . The DDFT is clearly not a theory for a colloidal fluid 
in which there is a temperature gradient. In order to incorporate such an effect in a 
microscopic theory, one would need to construct an external potential that couples to 
both the position and momentum degrees of freedom of the colloids. In reality, such 
thermal effects normally also influence the solvent particles. In the present description 
we have effectively assumed that the solvent acts as a heat bath, keeping the temperature 
of the colloids at a constant value - even in cases where the change in the external 
potential is such that one would find an increase in the temperature of the fluid were it 
a simple (molecular, non-colloidal) fluid - for example under rapid compression. 

As we have already mentioned, F ea .[{p g (r)}], the equilibrium excess Helmholtz free 
energy functional for the fluid mixture is, in principle, an unknown quantity. There 
exist, depending on the fluid in question, a number of accurate approximate functionals, 
which give extremely good results for equilibrium fluid density profiles for a wide variety 
of external potentials (see, for example, Ref. |2H| and references therein). This means 
that the DDFT (|20p is a very appealing theory for the dynamics of an inhomogeneous 
colloidal fluid, since it builds directly upon EDFT, one of the most successful theories for 
the equilibrium properties of inhomogeneous fluids. As the presentation above shows, 
the DDFT is clearly not an exact theory. However, the fact that, in principle, the 
equilibrium profiles obtained from the DDFT are almost exact leads one to expect that 
the theory should be reliable - at least when the fluid is not too far from equilibrium. 

3. Application of the DDFT to the GCM 

We argued in the previous section that the DDFT (|2T)|) should be a good theory for the 
dynamics of a colloidal fluid when the fluid is near to equilibrium, provided we have an 
accurate approximation for the excess Helmholtz free energy functional F ex [{pi(r, t)}] 
for the fluid (we now use i, j, rather than q, q' to label the different species of colloids). 
In this section we shall demonstrate that in the case of a particular model fluid for 
which we do have an accurate Helmholtz free energy functional, the DDFT seems to be 




(22) 




(23) 



On combining this result with the continuity equation 




(24) 
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reliable even for situations when the fluid is quite far from equilibrium and phenomena 
such as phase separation and interfacial adsorption (wetting) are present. Throughout 
we assume that the Brownian level of description Q provides an accurate account of 
the underlying particle dynamics. 

Dzubiella and Likos ^T] applied the DDFT to a one-component model fluid in 
which the particles interact via a purely repulsive Gaussian potential, the Gaussian 
core model (GCM). Their choice of model fluid was motivated by the fact that a 
simple mean-field approximation for F ea; [{/9j(r)}] proves to be quite accurate for the 
equilibrium properties, and thus they were able to test whether the DDFT formulation 
(llfjjl of dynamics is accurate. We shall say more about this model fluid below. Their 
strategy was to consider the inhomogeneous fluid confined in either a spherical cavity or 
a slit. They used EDFT to calculate the equilibrium density profile corresponding to a 
particular external potential and considered cases when the external potential suddenly 
changed (i.e., a parameter in the external potential was either increased or decreased). 
Using the density profile from the EDFT as the starting density profile, they used the 
DDFT to determine how the density profile of the fluid evolved towards equilibrium. 
The reliability of their results was assessed by making comparison with BD simulation 
results - i.e. they numerically integrated Eq. (^J) a number of different times, for different 
realisations of the stochastic noise, and then averaged over all the different runs in order 
to obtain the ensemble average time evolution of the fluid density profile. Dzubiella and 
Likos found that the DDFT and the BD simulation results were in very good agreement 

We follow the same strategy here but for a binary mixture of GCM particles. 

The one component GCM does not exhibit fluid-fluid phase separation [3J[Tn|. It is 
therefore of interest to find out whether the DDFT proves to be reliable for the dynamics 
of inhomogeneous fluids when there is phase separation in the fluid and when related 
interfacial phenomena such as wetting of the cavity wall are present. The model fluid 
we consider is the binary GCM. This fluid does exhibit liquid-liquid phase separation 
|2U[ I21j. The GCM particles interact via the following purely repulsive pair potential: 



and no other higher body potentials. In Eq. ()25|) i, j = 1, 2 label the two different species 
of particles, e^- > is a parameter that determines the strength of the interaction and 
Rij denotes the range of the interaction potential. Note that this potential has no 
hard core; the centres of the particles can overlap completely. When one considers 
the effective potential between the centres of mass of polymers in a good solvent, one 
finds that the Gaussian potential (J23j) . with e^j ~ 2k B T and R it j ~ R g , the polymer 
radius of gyration, provides a good approximation [31 H2J ECU UH ESI EE] • When the fluid 
density becomes sufficiently high, each GCM particle interacts with a large number of 
neighbours, and it is established that the following mean-field excess Helmholtz free 
energy functional [3 EH 1201 HQ : 



u iA r ) = %i exp(-r 2 : / 'R^) 



(25) 




(26) 
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becomes rather accurate. This functional generates the RPA closure for the direct pair 
correlation functions (c.f. Eq. ((7j)): 

^-^WwT- 8 "^-^ (27) 

We solve the DDFT (J2(J|) for the binary GCM confined in spherically symmetric 
external potentials of the form: 

uUr) = E(r/TZ) 10 , (28) 

where r is the distance from the origin, E = lOkgT and the length-scale TZ is the same 
for both species of particles. This potential is of the same form as one of the external 
potentials considered in Ref. |TI] for a one component fluid of GCM particles. As the 
length parameter TZ is increased, the size of the cavity is increased. We consider cases 
where for times t < the fluid is at equilibrium confined in a cavity with potentials f!28|) . 
with TZ — TZ\. Then at t — the cavity potentials change to ones with TZ = TZ 2 7^ TZ\. 
Due to the spherical symmetry of the external potentials, the (ensemble average) fluid 
one body density profiles will also display spherical symmetry. In the binary GCM fluid 
we shall consider mixtures with two different sets of pair potential parameters. The first 
corresponds to a set giving bulk liquid-liquid phase separation and the second to a set 
giving microphase-separation [23 . We shall also assume throughout that the friction 
constants for the two different species of particles are equal, i.e., T^ 1 = T^ 1 = P -1 , to 
keep the problem as simple as possible.^ 



3.1. DDFT for a GCM fluid which exhibits bulk phase separation 

We consider a binary GCM fluid with pair potential parameters en = e 22 = 2/c#T, 
e 12 = 1.8877fc B T, R 22 = 0.665i? n and R 12 = 0.849i? n , which exhibits bulk fluid-fluid 
phase separation. The fact that R\ 2 = (1 + A)(i?n + R 22 )/2, with A > 0, ensures 
that the fluid exhibits positive non-additivity and it is this feature which induces phase 
separation. This choice of pair potential parameters was used in a number of previous 
studies by the author [2TJ 122 HHH EU 05] • The fluid phase-separates at high total 
densities p = p\ + p 2 , where p\ and p 2 are the bulk densities of the two species. Within 
the mean-field DFT defined by Eq. (|2T)j) the (lower) critical point is at pR\ x = 5.6 and 
concentration x = p 2 /p = 0.70 [2"T] . 

The first case we consider is that with Aq = 200 big particles of species 1 and 
N 2 = 100 particles of species 2, confined in a cavity with TZ = 4i?n for t < 0. At t = 
the cavity size is suddenly reduced to TZ = 3Rn. In Fig. [Uwe display the evolution of 
the fluid density profiles after this sudden compression of the cavity. We display the 
results from solving the DDFT (|2"U|) combined with the RPA functional 1)26)1 . as well 
as results from BD simulations j^SJ. In the BD simulations we numerically integrate 
Eq. © generalised to two different species of GCM particles. We typically perform 

If Another more realistic choice, following Stokes, would be to set iy 1 oc Ru. However, we do not 
believe making this alternative choice would affect any of our overall conclusions. 
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Figure 1. Density profiles Pi(r,t) (DDFT results: solid line for species 1, dashed 
line for species 2; symbols are BD results, (+) for species 1, (x) for species 2) for a 
fluid composed of Ni = 200 particles of species 1 and N 2 = 100 particles of species 
2, which is initially (t < 0) at equilibrium in an external potential of the form in Eq. 
(|28H with TZ = 4Rn. At t = the external potentials suddenly change to those with 
TZ = 3i?n. The profiles are plotted for various t* — kBTTR^t. The two species of 
particles are uniformly mixed in the cavity at t — 0, but due to the increase in density, 
the equilibrium profiles for the fluid in the cavity with TZ = 3i?n exhibits a degree of 
demixing. Note also that p\{r — 0,t) is a non-monotonic function of time. 




012345 012345 012345 

r/R„ r/R„ ;-/R„ 



Figure 2. This is the reverse case of that displayed in Fig.^ Initially (t < 0) the fluid 
is at equilibrium confined in external potentials of the form Eq. (|28|) with TZ = 3Rn. 
Then at t = the cavity potentials suddenly change to those with TZ — 4i?n. The 
profiles are plotted for various t* — ksTTR^t. 
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Figure 3. The second moment of the density profiles, m\{t), defined by Eq. J2HJl. 
The left hand figure is for the cases corresponding to the density profiles displayed in 
Figs. H and [21 The right hand figure corresponds to the cases in Figs. 01 and The 
solid lines are the second moment for species 1 , m\(t), and the dashed lines m%(t), for 
species 2. In each case, the curves with a higher value of m l 2 (t — > oo) correspond to 
a final external potential with a larger value of TZ. Note that in both cases m\(t) is 
non-monotonic, whereas m\(t) is a monotonic function of time. 



1500 different runs, each with different realisations of the stochastic noise term, and 
then average over all the different runs in order to obtain the ensemble average time 
evolution of the fluid density profiles. In order to generate each starting configuration, 
we allow the fluid to equilibrate for a reduced time t* = ksTTR^t = 10 with a fixed 
parameter TZ = TZ\(= 4Rn) in the external potential, before changing this parameter 
to TZ — TZ 2 (= 3i?u) and determining the relaxation of the fluid. In Fig. [J the density 
profiles at t — correspond to those of the equilibrium fluid in a cavity with TZ = 4i?n + . 
When the external potential parameter TZ = 4Rn, the two species of GCM particles are 
mixed within the cavity, although because the walls of the cavity favour species 1, there 
is a higher density of species 1 around the outside of the cavity. This preference of the 
cavity wall for species 1 is due to the fact that the cavity wall potentials for both species 
of particles decay into the fluid with the same decay length. This results in an effective 
attraction between the wall and species 1 [22J. For small t > 0, the initial 'compression' 
of the fluid results in the growth of a sharp peak in the densities of both species around 
r ~ TZ2 = 3i?n- This causes a density 'wave' to travel through the fluid into the centre 
of the cavity. The fluid reaches equilibrium at t* = t/r B ~ 1.3, where tb = (3/YE^ l is 
the Brownian time scale. In Fig. ^ we display the fluid density profiles for the times 
t* = 0, 0.025, 0.1, 0.3, 0.8, and 2. The profile for t* = 2 effectively corresponds to the 
equilibrium profile for the external potential ()28|) with TZ = 3Ru. With this value of TZ 
the fluid exhibits a degree of phase separation in the cavity due to the fact that the total 

+ Note that within the grand canonical EDFT the mean number of particles (iVj) of species i in the 
cavity is constrained. In practice we set the chemical potentials fii such that the the mean number of 
particles in the cavity, (iVj) = Jdrpj(r) = iV», the number of particles in the BD simulation. 
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Figure 4. Density profiles pi(r,t) (solid line for species 1, dashed line for species 2) 
for a fluid composed of N\ = 25000 particles of species 1 and N 2 = 5000 particles of 
species 2, which is initially (t < 0) at equilibrium confined in external potentials of 
the form in Eq. H28|) with 1Z — 9i?n. Then at t = the external potentials suddenly 
change to those with 1Z — 13i?n. The profiles are plotted for various t* = ksTrR^t. 
Initially, the fluid is separated into two phases, one rich in species 1 around the outside, 
'wetting' the wall of the cavity, and the other phase in the centre of the cavity, rich in 
species 2. Notice the slow diffusion of particles through the fluid-fluid interface in the 
later stages of the equilibration. 



density of the fluid in the cavity has increased. The particles of species 1 are mostly to 
be found adsorbed around the wall of the cavity and species 2 towards the centre of the 
cavity. In Fig. |21 we display the density profiles for the reverse situation to the case in 
Fig. [TJ The external potentials (|2~%j) are initially (t < 0) those with 7Z = 3Rn. Then, at 
t = 0, the potentials are changed to those with 71 = 4i?n. In this case the fluid reaches 
equilibrium when t* ~ 1.7. The reason that the 'compression' case (Fig. is able to 
reach equilibrium faster than the relaxation case (Fig. EJ) is that in the compression 
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Figure 5. This is the reverse case of that displayed in Fig.0J Initially [t < 0) the fluid 
is at equilibrium confined in external potentials of the form in Eq. i|28[) with 1Z = 13i?n. 
Then at t = the external potentials suddenly change to those with 1Z = 9Rn. The 
profiles are plotted for various t* — kgTTR^t. Notice the concentration 'wave' (dip 
in the profile of species 1 and a peak in profile for species 2) that travels towards the 
centre of the cavity which allows the fluid to reach equilibrium slightly faster than in 
the opposite (relaxation) process - see Fig. ^ 



case the 'shock wave' that travels inward, mixes up the fluid, which allows it to reach 
equilibrium faster than in the relaxation case. In Fig. El we display the second moment 
of the density profile, 

m* 2 (t) = J drr 2 Pi (r,t). (29) 

Dzubiella and Likos JT] found that for a one component GCM fluid this is a monotonic 
function of time. For the relaxation (cavity expansion) case they found that m2(t) 
can be very accurately approximated by an exponentially decaying function. For the 
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compression case, m 2 (t) is also monotonic, and can be accurately parameterised by a 
function composed of two exponentials. They also found that the one component GCM 
fluid reaches equilibrium faster in the compression case, where TZ is decreased at time 
t = 0, rather than the relaxation case. For the binary GCM fluid, m\{t) is also a 
monotonic function of time (dashed lines), but interestingly, m\(i) is not (solid lines) - 
see Fig. El This is because particles of species 1 are adsorbed on the wall of the cavity 
and when the parameter TZ in the external potentials is suddenly decreased, particles 
of species 1 are forced more strongly towards the centre of the cavity than the species 
2 particles. This causes an increase in the density of species 1 at the centre of the 
cavity. However, in the final stages of the fluid equilibration, the density of species 
1 at the centre of the cavity decreases again. In other words, pi(r = 0,t) is a non- 
monotonic function of time. Similarly, when the cavity is suddenly increased in size, it 
is the particles of species 1 that 'feel' the space around the outside of the fluid that has 
suddenly appeared, rather than the species 2 particles, and so it is species 1 particles that 
move outwards to fill this space. However, in the final stages of the fluid equilibration 
process, there is a net flow of species 1 particles back towards the centre of the cavity. 
For these reasons, only m\(t) is a non- monotonic function of time. Dynamic processes 
such as these only occur when the fluid confined in the cavity is a binary fluid, where 
the wall of the cavity has a preference for one of the species. As can be seen from Figs. 
Q arid El there is remarkably good agreement between the density profiles obtained from 
the DDFT and from the BD simulations. This agreement gives us confidence concerning 
the reliability of the DDFT for cases where the numbers of particles are such that BD 
simulations become computationally too expensive. Note also that the RPA functional 
(f26j) becomes increasingly accurate as the GCM fluid density is increased 01113; this 
should further improve the accuracy of the DDFT results. 

In Figs. El and El we display the DDFT results for a similar situation as in Figs. 
Q and El except in this case the cavity is larger and the number of confined particles 
is higher: N\ = 25000 particles of species 1 and N 2 = 5000 particles of species 2 (the 
average densities are also higher). In Fig. Elwe display the results for the case when 
the fluid, for t < 0, is at equilibrium in a cavity with potentials given by Eq. (|2*H|) with 
TZ = 9Rx\. Then at t = the potentials change suddenly to those with TZ = 13i?n. 
We plot the density profiles for t* = 0, 0.02, 0.04, 0.06, 0.1, 0.5, 1, 5, 10, 20, 30 and 
oo. (Here and elsewhere, the t = oo density profiles are those obtained using EDFT 
for cavity potentials with TZ = TZ 2 . In the present case TZ 2 = 13i?n.) We see that at 
t = the fluid is strongly phase separated in the cavity, with the phase rich in species 
1 'wetting' the wall of the cavity, and the phase rich in species 2 at the centre of the 
cavity. The final equilibrium configuration, t — > oo, when the cavity radius TZ = 13 Rn, 
is that where the two species of particles are mixed together in the cavity, although the 
preference of the cavity wall for species 1 ensures that the density of species 1 is still 
higher around the outside of the cavity. In order to reach this equilibrium configuration 
the fluid first exhibits a 'quick' flow of the phase rich in species 1 to fill the space created 
around the outside of the fluid by the cavity expansion. There is then a second 'slow' 
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process whereby the two demixed phases mix, i.e. the diffusion of particles across the 
fluid-fluid interface is a slow process. 

For the reverse (cavity compression) situation, the results are displayed in Fig. El 
For t < the fluid is at equilibrium in a cavity with potentials given by Eq. (|28|) with 
1Z = 13i?n . Then at t — the potentials change suddenly to those with 1Z = 9Rn. 
We plot the density profiles for t* = 0, 0.02, 0.1, 0.5, 1, 3, 5, 10, 15, 20, 35 and oo. In 
this case the fluid reaches equilibrium slightly faster than in the case where the cavity 
size is increased. This is because in the cavity compression case, the sudden decrease of 
1Z at t — sends a particle concentration 'wave' into the centre of the cavity (see Fig. 
EJ) which allows the fluid to reach equilibrium faster. This 'wave' forms as a dip in the 
density profile of species 1 and a peak in profile for species 2 couple together and move 
towards the centre of the cavity. The amplitude of the wave increases with proximity 
to the centre of the cavity. As can be seen in the right hand figure of Fig. El for both 
these cases the second moments of the density profiles, given by Eq. (}2*9*|) . show similar 
behaviour to the cases with fewer particles in Figs, d and El see Fig. El m\{t) is a 
non-monotonic function of time, whereas rn\(t) is a monotonic function of time. As for 
the cases in Figs. [Hand El this is due to the fact that the species 1 particles are adsorbed 
on the cavity wall, rather than the species 2 particles. 

3.2. DDFT for a G CM fluid which exhibits microphase separation 

We now consider a binary GCM fluid with pair potential parameters the same as in 
the previous subsection except that i? i2 = 0.6i?n, i.e. R\ 2 = (1 + A)(i? n + R 2 2)/2, 
with A < 0; the fluid exhibits negative non-additivity. This fluid does not exhibit bulk 
fluid-fluid phase separation. Due to the negative non-additivity there is a propensity to 
ordering in the fluid in which particles of species 1 have as nearest neighbours particles 
of species 2, and vice versa [23! ■ We refer to this phenomenon as '1-2 ordering'. At 
high total densities EDFT predicts that this GCM fluid freezes into a crystal in which 
the particles are highly delocalized, with Lindemann parameters as high as 90% near 
melting [231- When the fluid is confined in a spherical cavity and the fluid density 
is high enough, the fluid forms an 'onion' structure of alternating layers, one particle 
thick, of the two different species [23! ( see a ^ so Figs. Ifij [7[ 191 and HP} ■ Similarly, when 
the high density mixture is confined in a planar slit the density profiles show that fluid 
forms alternating layers of the two different species parallel to the walls of the slit [23J. 
At lower densities the two different species of particles are mixed. This microphase 
separation is associated with the fact that the bulk fluid exhibits an instability with 
respect to periodic density modulations - a 'A-instability'. For more details regarding 
the origin of this instability see Ref. [23J and references therein. Here our interest is 
limited to the question: Can the DDFT describe the formation of 'onion' structures in 
a spherical cavity with potentials given by Eq. (j2*H|) when 1Z is reduced, and can the 
DDFT describe the onion 'melting' when TZ is increased? 

We follow a strategy similar to that taken in the previous subsection and make 
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Figure 6. Density profiles pi(r, t) (solid line is DDFT results for species 1, dashed line 
for species 2; symbols are BD results, (+) for species 1, (x) for species 2.) for a fluid 
composed of iVi = 200 particles of species 1 and N 2 = 150 particles of species 2, which 
is initially (t < 0) at equilibrium confined in external potentials of the form in Eq. 
(12811 with 1Z = 4i?n. Then at t = 0, the external potentials suddenly change to those 
with 1Z = 3Rn. The profiles are plotted for various t* = ksTTR^t. This model fluid 
exhibits microphase-separation. The final t* = 1 configuration is an 'onion' structure. 




Figure 7. This is the reverse case of that displayed in Fig. El Initially (t < 0) the fluid 
is at equilibrium confined in external potentials of the form in Eq. I|28|) with 1Z = 3i?n. 
Then at t = 0, the external potentials suddenly change to those with 1Z = iRn. The 
profiles are plotted for various t* = ksTYR^t. The initial t — configuration is an 
'onion' structure. 
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a comparison between the results from DDFT and BD simulations for cases with a 
limited number of particles, in order to demonstrate the reliability of the DDFT for 
this particular GCM fluid and then we apply the DDFT to cases where the number 
of particles is too high for simulations to be a realistic tool for studying the dynamics. 
In Fig. E] we display the density profiles for a fluid composed of N\ = 200 particles of 
species 1 and N 2 = 150 particles of species 2, which is initially (t < 0) at equilibrium 
confined in external potentials of the form in Eq. ()28|) with TZ = 4i?n. Then at t — 
the external potentials are suddenly changed to those with TZ = 3Rn. We display the 
density profiles for t* = 0, 0.025, 0.05, 0.1, 0.5 and 1. In the initial t = density profiles 
we see that there is little sign of the fluid exhibiting 1-2 ordering. As was the case for 
the fluid described in the previous subsection, the wall potentials result in an effective 
attraction between the wall of the cavity and particles of species 1. Therefore there is an 
increased density of species 1 around the wall of the cavity. Following the compression in 
the cavity at t — 0, there is an initial in-flow of fluid towards the centre of the cavity and 
an increase in the density of the fluid near to the cavity wall. Following this initial stage, 
the fluid relaxes into a configuration which exhibits a pronounced degree of 1-2 ordering 
- see the t* = 1 figure in Fig. El In Fig. Owe display the results for the reverse situation, 
i.e. when the fluid is initially (t < 0) at equilibrium confined in external potentials of 
the form in Eq. (|2*B|) with TZ = 3i?n, then at t = 0, the external potentials suddenly 
change to those with TZ = 4Rn. We display the profiles for t* = 0, 0.025, 0.05, 0.1, 
0.5 and 1. In both cases the agreement between the DDFT and the BD simulations is 
remarkably good. Fig. |H1 displays the results for the second moment of the density profile, 
defined by Eq. (j2HJ). As in the previous subsection, in the present case the moments 
m\(i) for species 1 are non-monotonic functions of time, whereas the moments m\(t) 
are monotonic functions of time. The fluid which exhibits 1-2 ordering is able to reach 
equilibrium faster than the fluid considered in the previous subsection which exhibits 
bulk phase separation. This is because the average distance the particles must diffuse to 
be arranged in a microphase-separated distribution is shorter than in the case where the 
final equilibrium configuration is that exhibiting 'bulk' phase-separation. In the latter 
case the particles of species 1 must diffuse to the outside of the cavity, whilst particles 
of species 2 must diffuse to the centre of the cavity in order to reach equilibrium. 

In Figs. IH1 and El we display the DDFT results for a similar situation, except now 
the cavity is larger and the number of confined particles is also much bigger: Ni = 16000 
particles of species 1 and N 2 = 15000 particles of species 2 (the average densities are 
also higher). In Fig. El we display the results for the case when the fluid, for t < 0, is at 
equilibrium in a cavity with potentials given by Eq. (|2"%|) with TZ = 8Rn. Then at t = 
the potentials change suddenly to those with TZ = 10i?n. We plot the density profiles 
for t* = 0, 0.02, 0.06, 0.1, 0.2, 0.5, 0.7, 1 and oo. At t = the two different species 
of particles are strongly ordered into alternating layers of the two different species - an 
'onion' structure. The final equilibrium configuration, t —>■ oo, is that where the two 
species of particles are uniformly mixed in the cavity, although, due to the preference 
of the cavity wall for species 1, the density of species 1 is higher in the outer region of 
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Figure 8. The second moment of the density profiles, m\{€), defined by Eq. Ij29(l . for 
a GCM fluid exhibiting microphase-separation, with density profiles displayed in Figs. 
El and H The solid lines are the second moment for species 1, m\{t), and the dashed 
lines m|(i), for species 2. In each case, the curves with a higher value of m\(t — > oo) 
correspond to a final external potential with a larger value of 1Z. Note that in both 
cases m\(t) is non-monotonic, whereas m\(t) is a monotonic function of time, as was 
the case in Fig. for a GCM fluid that exhibits bulk phase separation. 

the cavity. The equilibration process is almost an order of magnitude shorter in time 
than in either of the cases in Figs. 0] or for the GCM fluid exhibiting bulk liquid-liquid 
phase separation, even though the number of particles involved in these cases is similar. 
This is because the particles must diffuse a distance ~ Rn, the particle size, in order 
to go from a state with 1-2 ordering to a mixed state, whereas for the fluid exhibiting 
'bulk' phase separation the particles must diffuse a much larger distance ~ 7Z, the radius 
of the cavity. For the reverse (cavity compression) situation, the results are displayed 
in Fig. El Initially, for t < 0, the fluid is at equilibrium in a cavity with potentials 
given by Eq. (|2~%j) with 7Z = 10 Rn. Then at t — the potentials change suddenly to 
those with 71 = 8R U . We plot the density profiles for t* = 0, 0.01, 0.02, 0.06, 0.1, 
0.2, 0.5, 0.7 and oo. As with the cavity expansion case in Fig. El this process is almost 
an order of magnitude quicker in time than in either of the cases in Figs. 0] or El The 
compression causes the total fluid density to first increase around the wall of the cavity, 
before increasing near the centre of the cavity. This results in the 'onion' layers forming 
first in the outer region and then developing inwards towards the centre of the cavity as 
the fluid equilibrates. 

4. Liouville, Kramers and Smoluchowski equations 

We now turn to the question of what approximations are involved in the present 
description of the fluid dynamics. Under what circumstances can the equations of 
motion for a solution of colloids suspended in a solvent of much smaller particles be 
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Figure 9. Density profiles pi(r,t) (solid line for species 1, clashed line for species 2) 
for a fluid composed of Ni — 16000 particles of species 1 and N2 = 15000 particles of 
species 2, which is initially (t < 0) is at equilibrium confined in external potentials of 
the form in Eq. I)28|l with 1Z = 8Rn. Then at t = the external potentials are suddenly 
changed to those with 1Z = 10i?n. The profiles are plotted for various t* = huTTR^. 
This model fluid exhibits microphase-separation. The initial t = configuration is an 
'onion' structure. 

approximated by the stochastic equations of motion (jHJ)? Or, equivalently, in which 
situations can the time evolution of the probability density function for the full system 
of colloid and solvent particles be approximated by the Smoluchowski equation (JTTJ)? 
In Refs. [37J EE] these questions are addressed for the case of a fluid in which one 
of the particles is much larger than the rest and when the solvent friction coefficient 
r _1 is sufficiently large, that one can argue that the time evolution of the probability 
distribution function for the single big (colloid) particle is governed by the single 
particle Smoluchowski equation. Here we present an argument which is a simplified 
generalisation (for N big colloid particles) of this argument. The derivation does not 
contain new ideas. However, we do make connections between older, well known, results 
concerning the dynamics of colloidal fluids and more recent developments in the theory 
of solvent mediated effective potentials between colloids in solution [3 . The following 
therefore provides much insight to the physics incorporated in the Smoluchowski 
equation and therefore also in the DDFT, and applies generally to colloidal fluids. 

We consider a fluid consisting of a single species of N colloid particles of mass m, 
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Figure 10. This is the reverse case of that displayed in Fig. Initially (t < 0) 
the fluid is at equilibrium confined in external potentials of the form in Eq. H28JI with 
1Z = lCLRn. Then at t = 0, the external potentials suddenly change to those with 
1Z = 8Ri%. The density profiles (solid line for species 1, dashed line for species 2) are 
plotted for various t* = ksTTR^t. 



suspended in a single component solvent composed of S solvent particles, of mass M. 
We denote the coordinates of the i th colloid by and momentum pj. The set of colloid 
position coordinates we denote by = {ri,...,rjv} and similarly = {pi, Pjv}- 
For the solvent particles we denote the location in phase space of the v th solvent particle 
by (Rj,,Pj,). Just as for the colloids, we denote the sets R 5 = {Ri,...,Rs} and 
P s = {Pi, Ps}- The Hamiltonian for this system is: 

i=l u=l 

+ V N (r N , t) + ^(r^, R s ) + V S (R S , t), (30) 

where the first two terms on the right hand side are the colloid and solvent kinetic 
energy contributions to the Hamiltonian. VN{r N ,t) is the colloid potential energy: 

N 1 N 1 N 

V N (r N ,t)=J2 v ext(ri,t) + -J2J2 V2 ^ r ^ + Q 5ZX^ 3 ( r *' r i' Ffc ) + - ( 31 ) 

i=l j^i i=l kjtjjti j^zi i=l 

which is assumed to be made up of a one-body term (external potential v ext (ri,t) on 
each colloid particle) a two-body term (u 2 (r i ,r :) ) is the pair potential), a three body 
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term f 3 (rj, r^, r^), etc. Note that Vjy(r , t) describes the direct or bare colloid-colloid 
interaction potential; it does not involve effective solvent mediated interactions and is 
therefore not the same as Uj^{r , t) in Eq. (JHJ) - we will further clarify this issue later 
(see Eq. (HHJ)) *. Similarly, Vs^R 5 , t) is the solvent potential energy and Vjv,s( rAr j R S ) 
is the potential energy arising from interactions between the colloid and the solvent 
particles. 

We can define a phase space probability density function f( N+s ^(r N , p^, R s , P s , t) 
for the fluid, and its time evolution is governed by the (exact) Liouville equation |29j : 

g^g+g i_A gff+g A gg+g _i_^ p cy (JV+s) 

dt + m Z^ Pi ' dTl + Z- Xi ' d Pi v dR v 

+ £v^t^,^ + £z,^0, (32 ) 

j/=l y=l i=l 

where 

Xl -"9r7' Y "--dR u > Zj -"^T' (33) 
are forces on the particles. The Liouville equation is a statement of the continuity of 
f{N+s) j n phagg space over time. 

Since the solvent particles are much smaller than the colloids, they will equilibrate 
on a time scale t s much smaller than the time scale r c on which the colloids equilibrate, 
i.e. r c ^> t s . We are interested in phenomena that occur on time scales ~ r c , so we can 
assume that effectively the solvent particles equilibrate instantaneously, and when we 
integrate over the solvent degrees of freedom in Eq. (J32j) we obtain an equation similar to 
the Liouville equation for a one component fluid of N particles, but with an additional 
'solvent' term due to interactions between the colloid and solvent particles: 

dfW 1_A ^ x dfW _ fdf^\ 

fit m ' P * fir ^— — ' ' fin I fit I 



(34) 



solvent 



dt m ^— ' dri ^— ' dpi V dt J 

i=l 1=1 s ' 

where the colloid reduced probability density function is 

fM(r N ,p N ,t) = J dV s J dR s f N+s \v N ,p N ,R s ,P s ,t) (35) 
and formally the 'solvent' term is 

The term / dP s f dR s Z i f ( - N+s ^ in is proportional to the average force exerted on 
colloid i by the solvent particles. On this 'coarse-grained' time scale ^> r s , the leading 
order contributions to this force are a one body (Stokes) drag force on each colloid, 

* For generality we assume Vn(t n ,t) contains higher body terms, although it is generally assumed 
that the bare interactions between particles are pairwise with no higher body terms. In principle, the 
higher-body terms arise from integrating over quantal degrees of freedom in obtaining Vjv(r , i). 
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— 7Pj, where 7 is a friction coefficient, and a force term, Xj, due to solvent mediated 
interactions between the colloids. This force is 

x, = (37) 

where <&n(t n ) is the effective solvent mediated potential between the colloid particles. 
For an equilibrium fluid $>n(t n ) = —k B T 'In Q s(r 1 ) , where Qs(r N ) is a partial partition 
function [3]: 



-35 



f 



{V N ,s(r N ,R S ) + V s (R s ,t)) 



(38) 



k B T 

and A5 is the de-Broglie wavelength for the solvent particles. For the equilibrium fluid 
$jv(r ) can, in principle, be calculated. In general 

N N N 

$ w (r Af ) = V0 o + ^0i(r i ) + -^^0 2 (r l ,r J ) + - £ £ £> 3 (ri, r.,-, r fc ) + ...,(39) 

i=i t=i fc^jV* j¥< i=1 

where K is the volume of the system, O is a zero-body potential j^j, 0i(i"j) is the solvent 
mediated one-body potential (for example with the fluid container walls), 02( r i, r j) is 
a two body pair potential, </> 3 (rj, r.,-, r fc ) is a three body potential, and so on. These 
potentials are generally density dependent. We shall assume that the non-equilibrium 
solvent mediated potential is the same as the equilibrium solvent mediated potential 
(|39|). This approximation should be reliable on time scales 3> t 8 , since on these time 
scales the solvent particles are effectively at equilibrium. Thus the leading order terms 
in a Taylor expansion of the force term in Eq. (j3T?j) can be expressed as: 

where 9 is a mobility coefficient, which, in principle, is a function of the colloidal phase 
space coordinates. However, we shall assume 7 and 9 to be constants. If we now 
substitute Eq. (j4*Uj) into Eq. retaining only the two leading order terms, then we 

obtain: 

d f(N) l " df (N) N {N) 

-^- + -22^-^— + z2 x *+ x * 

at m OYi ^— ' api 

i=l i=l 

^±±. fm + e ±^. (41) 

i=l ^ l i=l ^* 

We recognise this as the Kramers equation - i.e. a generalised Fokker-Planck equation 
[Sj. At this point we can also make the connection between the constants 7 and 9 via 
the fluctuation dissipation theorem, and we find that 9 = inksT^. In other words, 
when 9 = rnksTj, the equilibrium colloid reduced probability distribution function is 
correctly given by 



AN) 

Jo « ex P 



/ 1 N 




(42) 
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where Vn(t n ), the term describing the direct interactions between the colloids, is given 
by Eq. (}3~Tj) with a time independent external potential. 

We now clarify the difference between the potentials UN(r N ,t) in Eq. (JBJ) and 
V]y(r N ,t) in Eq. (f3*T|). We see that the total colloid potential UN(r N ,t), is given by 
(see Eqs. (13T)) and Pty ): 

U N (r N ,t) = V N (r N ,t) + <l> N (r N ). (43) 

i.e., we see explicitly here that UN(r N ,t) is the sum of the direct colloid interaction 
potential (j3~Tj) and the solvent mediated potential (|3Ti|). In Eq. (JHJ) we can therefore 
identify the total colloid effective one body potential u ext (ri) = v ext (ri) + 0i(i"j), as the 
sum of a direct potential and a solvent mediated potential. Similarly, the total colloid 
effective pair potential u 2 (ri, Tj) = v 2 (ri, Tj) + 2 ( r i, T j), the total colloid effective three 
body potential 1*3(1^, r^, r^) = ^(rj, Tj, r^) + 03 (r^, r 3 -, r^), and the higher body effective 
potentials are a sum of a direct contribution and a solvent mediated contribution. 

The Kramers equation (|41jl can also be obtained as the (generalised) Fokker-Planck 
equation for the time evolution of the probability density function for N colloid particles 
with the following stochastic equations of motion [H] : 
dr.^ _ pi 
dt m' 

It = ~ m " wy^^ + $7v(r7V)) + Gi(t) ' (44) 

where G«(t) = (£f (£), £f (£), £f(i)) is a white noise term with correlations given by 
Eq. (jl(Jj) . In a stochastic treatment of the colloidal fluid one considers the above 
equations of motion, (|44p. solely for the colloids and one incorporates the effect of the 
solvent via the stochastic noise term, the friction coefficient 7 and the solvent mediated 
potential <&n(t n ). In the treatment above, we arrived at the Kramers equation (14 lj) by 
approximating the Liouville equation for the full system of solvent and colloid particles. 

In deriving the Kramers equation (|41|) it is clear that we made a very large reduction 
in number of degrees of freedom in the description of the fluid, since we have integrated 
over the solvent degrees of freedom. However, there are still 6iV degrees of freedom 
in (|41|) and further reductions are required. We now focus on the colloid probability 
density function 

pW(rV) = / dp N fW(r N ,p N ,t). (45) 

In order to obtain an equation for the time evolution of P( N \r N ,t) we have to make 
further approximations. The following is a generalisation of the "quick and dirty" 
approach in Ref. [2] (see also Refs. jSJEn])- On integrating with respect to the colloid 
momentum degrees of freedom in Eq. 1)41)1 . we obtain the following continuity equation: 

+ E^- J ' = ' ( 46 ) 



dt ^ <9r, 

i=l 



where the current 



J 4 = I dp N -pJ^(r N ,p N ,t), (47) 
m 
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and we have used the fact that the integrals J dp^df^ /dpi), J dpi(d 2 / dpf) and 
J dpi(d(pif^ N ')/dpi) are all equal to zero. 

Also, if we multiply Eq. (|41jl through by pk/m, the velocity of the k th colloid, and 
then integrate over all the colloid momentum degrees of freedom, we obtain: 

y fd P NPk ( Pi 9fiN) ) i v [dp NPk (jib i xa df{N) ) 

J m \m dvi ) J m\ l 1 dp { J 



dt 



= 7E / dp N ^-pJ {N) -imk B Tj2 [ *P N -^T- (48) 
j^J mdpi j^J m dpf 

If we now make a 'local momentum equilibrium' approximation [2], which sets terms such 
as / dp iV Pi, a Pj,bf = fnksTP^SijSafi, where p^ a is the a-component (a, b = x,y,z) 
of the momentum of the i th particle, then we find that Eq. (|48|) reduces to 

dt m dr k m 
If the friction constant 7 is sufficiently large then we can neglect the first term in (@HJ), 
as it will be negligible compared to the friction term — 7 J*., on the Brownian time scales 
r B , and we obtain: 

j^-r^-r^pw, (50) 

OTfc or k 

where T = 1/7717. Substituting (|5Uj) into Eq. (j4*Uj) . we obtain the Smoluchowski equation 

dm- 

The Smoluchowski equation is generally presented from the stochastic viewpoint, 
starting from the Langevin equations of motion However, as shown above, one 
can argue for its use as an approximation to the exact Liouville equations for the time 
evolution of the probability density function for colloidal fluids. The above derivation 
demonstrates the physical conditions under which the time evolution of the probability 
density function of a fluid of colloidal particles is described by the Smoluchowski 
equation. By going to the Smoluchowski equation (with an effective potentials between 
the colloids), the description of the fluid is reduced to one based on only the position 
coordinates of the colloids, rather than the full set of phase space coordinates for the 
colloid and solvent particles, which is, of course, a significant simplification. 



5. Discussion and conclusions 

In Sec. |21 we introduced the DDFT of Marconi and Tarazona |H1 d for colloidal fluids, 
following a recent derivation in Ref. [H] of the DDFT starting from the Smoluchowski 
equation. We emphasise that this is an approximate theory. In Sec. H3 we applied the 
DDFT to the specific case of a binary fluid of GCM particles confined in a spherical 
cavity of variable size. We find that the DDFT is able to describe accurately the time 
evolution of the highly structured fluid one-body density profiles. We find a variety of 
collective dynamic processes that are accounted for by the DDFT. For example, the 
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case in Fig. where we see a particle concentration 'wave' travelling through the fluid 
into the centre of the cavity that allows the fluid to equilibrate into the phase-separated 
configuration slightly faster than one would expect, bearing in mind the results displayed 
in Fig. which showed that the fluid-fluid interface persists for a long time after the 
cavity radius 1Z is increased. Such phenomena are very difficult to simulate because of 
the large number of particles involved. It is in situations such as these where the DDFT 
should provide a useful tool to analyse the fluid dynamics. 

Our strategy is to first consider cases of only a few hundred particles where BD 
simulations can be used to assess the accuracy of the DDFT results. In these test 
cases with fewer particles, where the fluid exhibits a limited degree of phase separation 
and wetting behaviour, we find excellent agreement between the DDFT and the BD 
simulation results (see for example Figs. dandEJ). This gives us confidence to trust 
the reliability of the DDFT when applied to cases with many more particles, where 
comparison with simulation becomes computationally too expensive. We predict that 
were one to perform a BD simulation for (say) the cases presented in Figs. 0] and El 
one would find just as good agreement between the DDFT results and BD simulation 
results. In fact, one should argue that since the total fluid densities are higher in these 
situations, and the fact that the RPA functional ()26|) becomes more accurate at higher 
densities |H1I^], the results from the DDFT should be even more reliable in these cases. 

In Sec. EJ we derive the Smoluchowski equation, starting from the exact Liouville 
equations by making a series of approximations in order to go from a description of the 
fluid in terms of the full set of colloid and solvent phase space coordinates to a (far 
more manageable) description solely in terms of the colloid position coordinates (the 
Smoluchowski equation (fTTjl). Such integration over degrees of freedom is a process 
that is required in order to render a practical statistical description of any condensed 
matter system. We believe our derivation is useful because it highlights some of the 
approximations involved in the Smoluchowski equation for colloids and therefore also in 
the DDFT and thus sheds light on the status of the DDFT. 

The Smoluchowski equation is usually presented from the stochastic viewpoint, 
since it is the Fokker-Planck equation for the Langevin equations 0. This 
connection makes it clear that the Smoluchowski equation and thus the DDFT neglect 
hydrodynamic effects. To include hydrodynamic effects one would have to treat 
the solvent at a level beyond that taken here, where the effects of the solvent are 
incorporated in the DDFT via the effective solvent mediated potential between the 
colloids and in the friction coefficient One can build in hydrodynamic effects in the 
Smoluchowski equation by replacing V with a matrix Ty, describing the hydrodynamic 
coupling between colloids i and j jHJ EOJ |JT] • However, this makes the reduction of the 
Smoluchowski equation to a DDFT much less straightforward. 

We conclude that the DDFT, Eq. pOjl . should provide a good theory for the 
dynamics of the one body density of a colloidal fluid, even when there is phase 
separation and wetting of interfaces, provided that i) the friction constant characterising 
the solvent, is large enough and ii) there exists an accurate approximation for 
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the excess Helmholtz free energy functional T ex \p\. Many colloidal fluids can be 
accurately modelled as effective hard-sphere fluids and therefore we expect that the 
DDFT combined with fundamental measure free energy functionals |i2~l for hard 
spheres to provide a reliable theory for the dynamics of the one body density profile of 
a (hard-sphere) colloidal fluid. 
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